Almeida-Thouless transition below six dimensions 
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The existence of an Almeida-Thouless (AT) instability surface below the upper critical dimension 
6 is demonstrated in the generic replica symmetric field theory. Renormalization fiows from around 
the zero-field fixed point are investigated. By introducing the temperature and magnetic field 
dependence of the bare parameters, the fate of the AT line can be followed from mean field (d = oo) 
down to d = 6 — e. 
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Notwithstanding the relative simphcity of the relevant 
model postulated by Edwards and Anderson [l| , the Ising 
spin glass problem has resisted a thorough understand- 
ing for decades. Severe frustration makes numerical sim- 
ulations extremely hard and computer-time consuming, 
whereas analytical methods must handle the inhomogen- 
ities caused by the quenched disorder. The model was 
later extended and studied on the fully connected lattice 
by Sherrington and Kirkpatrick (SK) Q], the character- 
ization of the spin glass phase by the solution of Parisi 
(see Ref. Q for a list of references) is now unanimously 
accepted as the true mean field theory. This mean field 
spin glass proved to be very complicated, its eqilibrium 
state breaking up to ultrametrically organized ergodic 
components, commonly called pure states. This complex 
phase space structure survives in an external magnetic 
field up to a phase boundary called the Almeida-Thouless 
line. Approaching this line from the paramagnetic side 
an instability develops: using a replicated picture [ij, the 
diverging spin glass susceptibility signals the break down 
of the replica symmetric phase [J| , and replica symmetry 
breaking develops. 

An alternative theory — the so called droplet picture 



— emerged, however, and continues questioning the rel- 
evance of mean field ideas in finite dimensional systems 
0. In this theory the glassy phase is much simpler, and 
is limited to zero field: a convincing conclusion about the 
existence or lack of an AT line may resolve a decades long 
debate about the structure of the spin glass phase in the 
physical dimensions. Recent numerical simulations 0, 01 
in three dimensions essentially excluded the possibility 
of a transition in a field, whereas the four dimensional 
case remains somewhat ambiguous (see Q for references 
to earlier works). On the analytical side, we must men- 
tion the scaling considerations in Q and renormalization 
group (RG) calculations 13, U I , whereas a leading order 
field theoretical computation 1^ provided an AT line 
above 6 dimensions. This letter tries to dissolve the mis- 
belief that the AT line disappears below the upper critical 
dimension, by explicitly calculating it close to, but below 

Ising spin glass transition in an external magnetic 
field can be studied in the generic replica symmetric 
field theoretical model [l3| defined by the Lagrangean 
£ = £(2) + £i where 
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The n(n — l)/2 component fields are symmetric in the replica indices and (j)""' = 0; the spin glass limit requires 
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n ^ 0. The total number of spins N is included to ensure 
the correct thermodynamic limit, and momentum con- 
servation is understood in the primed summation. The 
zero-fi eld p aramagnetic phase corresponds to higher sym- 
metry |14{ , with all the bare parameters but mi and wi 
zero, it has a unique mass F, and below the upper crit- 
ical dimension = G the spin glass transition is gov- 
erned by the fixed point 
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■ = w* = e/(2 — n) and 
ml ~ —e/2, e ~ 6 — d. The mass is split by an external 
magnetic field into the three different components Fr, 
Fa and Fl, thus generating — while replica symmetry is 
still preserved — a kind of quadratic symmetry breaking. 
The crossover region is best studied by the introduction 
of the nonlinear scaling fields satisfying the exact 
renormalization flows (ji = Xi gi . [The Xi 's of the mass 
sector (i = 1,2,3) were computed from the renormaliza- 
tion flow equations, in leading order, in Ref. [l3|-] The 
RG equations provide a way to express the bare param- 
eters of the Lagrangean in terms of the scaling fields, 
hence the masses can be computed as functions of the 



Fr = ,9i-f2.g2+53 + 0(e). 

Fa = 31 - - 4)52 - (n - 3)53 + 0(e), 

Fl = 51 - 2(n - 2)52 + ^^^^^^.93 + 0{e). 



(3) 



The 0{e) terms neglected above have two contributions: 
the one-loop self-energy (which is computable at that 
order) and corrections to the bare masses expressed in 
terms of the giS — this, however, is not available in a 
leading order RG calculation. Coupling-like scaling fields 
5i's with z > 3 enter also at this 0(e) level. Eq. ([3]) 
can be derived by fixing the bare parameters such that 
((/)p'^) = 0; this condition determines 50 unambigously in 
terms of the other g^'s. Two critical surfaces can be found 
from Eqs. ([3]) in the low temperature (51 < 0) regime: 

• Fr = 0, Fa and Fl both positive — i.e. an 
Almeida-Thouless instability — for 52 < —51; 

• Fa = 0, Fl and Fr positive (n > 0) for 52 > —51. 

The common boundary of these two manifolds (which are 
two-dimensional now, but allowing for coupling-like scal- 
ing fields 5i, i > 3, they will have a complicated higher 
dimensional structure) for 52 = —51 is massive only in 
the longitudinal sector [ll| . 

We are now interested in the RG flows along the AT 
instability surface when starting in the crossover region. 
The first order RG equations were all presented in Ref. 
[l3|, their structure is best displayed by the following 
(temporary) redefinition of the couplings: Wi/^/e —f Wi, 
which are now, like the masses, order unity. With the 



scaling factor e'^', and t = el: 
dm. 



dl 

dwi 1 
~dr ~ 2 



2mi - eMi{mi,m2,m3;wi, . . . , ws), i = 1,2, 3; 

(4) 

m + Wi{mi,m2,m3;wi, . . . , Ws), i = 1, . . . ,8. 

(5) 



The Aii and Wi functions are quadratic and cubic, re- 
spectively, in the couplings. The most important feature 
of the RG equations above is that the flow parameter I 
in the mass sector, Eq. is much larger for e ^ 1 than 
t of the couplings, Eq. ([5]). Thus the masses renormal- 
ize in the background of the adiabatically slow couplings: 
the anomalous (A) and longitudinal (L) components, as 
they are 0(1) on the AT surface, blow up, whereas the 
replicon (R) one, 2mi = 0(e), evolves into its adiabatic 
fixed point determined by the initial values of the cou- 
plings and W2'^. While Wi~^ = ■wi{t = 0) = Wi, we 
must carefully follow the development of W2 in the tran- 
sient regime from W2{t — 0) = ^ to = W2{t) with 
e <C t ^ 1 for the following reason: For 1^1, i.e. t ^ e, 
wi and W2 decouple from the other bare parameters, and 
their flow can be put into the pair of equations: 



dwi 

~dr 

dr 



-wi +5„(r) wl 
-Kir) wl, 



where 5„ (r) and ft,„ (r) are cubic and quartic polynomials 
of r with coefficients which are simple polynomials of n, 
and r = 1x12 /wi. For the case rt = 0, these equations were 
derived and discussed in Ref. • We are now interested 
in the more generic case < ri ^ e, and observe that the 
qualitative behaviour of the renormalization flow changes 
drastically when the initial value of r(0) ~ Wj^/wJ^ 
passes through r\ — j^n -f O(n^), the unstable flxed 
point r\ being the solution of the equation hn{r) = 0. 
For r(0) > r^, we have runaway trajectories already no- 
ticed in [l^l with — > 00 and r ^ = 14.4 -|- 0(n); 
whereas for r(0) < rl, W2 immediately becomes negative, 
which is physically nonsense. 

To get r(0), we must integrate Eq. ([5]) for i = 2 in the 
transient regime from i = Otoe<ICi^l, thereby elim- 
inating nonrcplicon modes in the process of hardening 
anomalous and longitudinal masses. This is feasible us- 
ing Eq. (65) of Ref. ^ together with the table between 
the different sets of couplings in Eq. (49) of Ref. [l^l, 
resulting in r(0) < ior < n ^ 1 and starting close 
enough to the zero-fleld flxed point, whereas the spin 



^ More precisely, W2{t = 0) ^ e in the physically relevant part of 
the AT surface. 
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glass case is exceptional with the condition r(0) > 
always fulfilled. 

Runaway flows along critical surfaces have been asso- 
ciated with first order transitions in some common situ- 
ations with crossover phenomena [isj . and although this 
scenario cannot be ruled out completely for the spin glass 
either, we will argue that renormalization of the bare cou- 
plings on the AT surface towards their low-temperature 
limit may cause the runaway trajectories in this renor- 
malization scheme. To see this, we recall the derivation 
of the microscopic Lagrangean in Ref. [l3l |. and the ne- 
cessity to redefine the fields as c(j)p^ with c ^ T, 
to ensure the proper normalization of the kinetic term in 
C^^\ This will cause the couplings diverge even if they 
disappeared for T —^ otherwise. That kind of normal- 
ization was essential in the derivation of Eqs. ([4]) and (O, 
manifested in the flowing 77 exponents of the three differ- 
ent mass modes. As our approximate RG equations are 
valid only for Wi = 0(1), one probably needs to modify 
the RG scheme for detecting the proper zero-temperature 
behavior on the AT surface in this small e regime. This 
is, however, out of the scope of the present work. 

In the remaining part of this letter we want to locate 
the AT-line of the original Edwards- Anderson spin glass 
model on the AT-surface of the field theory above. For 
this reason, we must find out the dependence of the bare 
parameters in Eqs. ^ and ^ on temperature (T) and 
magnetic field (H). The criterium which is adopted here 
is that the tree approximation of the field theory (i.e. 
neglecting loops) be equivalent with the accepted mean 
field theory of the Ising spin glass, the SK model, whose 
replicated partition function has the form 0: 
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^ dqafi) and J^, the variation of 



the Gaussian distribution of the random Ising interac- 
tions, sets the energy scale. In the tree approximation 
fluctuations are omitted, which can be achieved by set- 
ting 0p£o = VN qap and zero for 4>p% in HI) and (g]), 
and comparing it with ([6]) and ([7]). Not forgetting that 
the bare parameters of the field theory are finally tuned 
by the transformation c/jp'^ — ^/N q S^Lq — *■ (j)"^^ — render- 
ing the one-point function to zero — , where q is the exact 
replica symmetric order parameter, they are expressed by 
T and H in the vicinity of the mean field critical point 



kT^^ = J as follows 
1 

wh 



\{H/kTf ~ (mic - T){wq) + i(n - 2){wqf + . . 



1 



uoi + U02+ 



mi = (mic - t) + (luq) + {H/kTY - ^ 9 

-{n - 1) uo3 + -^n{n - 1) U04 
TO2 = -{wq) - {H/kTf - ^ 9^ [(" - 3) uoi + W03] + • ■ ■ 



(8) 



^3 = ^q1^ ["01 + 2^04] + ■ • ■ 

ivi = w + 0{q,H^), and w, = 0{q,H'^) i = 2...8. 

Neglected terms above are higher orders in q and H^. 
The quartic couplings arc not included here, although 
their calculation is similarly straightforward, and they 
may be inportant above 8 dimensions, r > measures 
the distance from the critical temperature of the field 
theory (T^), whereas mic = - ^ {T^" ~ T^) gives the 
shift in the critical temperature, and is therefore one- 
loop order. The field theory is defined by r, (H/kT)^ 
and by the bare parameters of the symmetrical theory 
(zero magnetic field paramagnet): w (cubic coupling), 
"01, U02, uo3, (quartic couplings), . . .etc. (see [l9| for 
the classification of the quartic couplings) . To reproduce 
the SK model results in the tree approximation, we must 
put w = 1, uoi = 3, U02 = 2, uq3 ~ —6 and uoi = 0. 

The condition (^p'^) = provides us the equation of 
state, i.e. the order parameter q around T^, it is used here 
to eliminate r from our results, replacing it by q. The 
calculation of the one-loop contribution to the equation 
of state, and to the replicon mass is somewhat lengthy 
due to the complicated replica structure even in the case 
of replica symmetry. Nevertheless, it is still feasible by 
the methods of Ref. [l^. The result valid for d > 8, 
including the SK model by simply taking d = 00, can be 
put into the scaling form 



^ iwq)'^f^,{x,y), x 



[H/kTY 
[wqY 



and y ; 



[wq) 
(9) 



The scaling function Fr has the simple linear form 

T-ii{x,y) ^ ax + by + c, d > 6, (10) 

with a and b analytical down to 4 and 6 dimensions, 
respectively, and having their loop-expansions in terms 

a = 1 — 2w^ li and 

(11) 

6 = 1- 2ui^/6 + (-M10 + W30 + 4u4o)/4; d > 6. 

c however blows up at 8 dimensions, due to the infrared 
divergence developing in the first order contribution be- 
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hind the mean field term: 
2 

c = — —U2ow~^ ~16w^ Is + terms with /g and I4, d> 8. 

(12) 

As a resuh, scahng of the rephcon mass turns to the 
fohowing form when 6 < d < 8: 



{H/kTf 
with X = - — r-rj^— j- and y ; 



(13) 



{wqy/ 



2-3 ■ 



The scahng function preserves the form (jlOp with a and 5 
in (fTT|) . the constant c however becomes, instead of (fT2|l : 



-1610^ 



°° rf'^p 1 
(2^^ripT2)2 



, 6<d<8. (14) 



The zeros of the scaling function provide the AT tran- 
sitions, and two important cases can be studied for c? > 6: 

• = 0, i.e. a; = and y = —c/b. This case has 
been discussed in [l9| . 

• The spin glass limit 71 = 0, i.e. y = Q and x = x^) = 
—c/a. The AT line close to Tc in the two regimes 
is: 

UH/kTf = xo{wqf 8<d, 

[(i7/fcr)2 =a;o(wg)'^/2-i 5 < d < 8. ^ ^ 

From (HI]) and ([12]), with M20 = 2, the SK value xq = 4/3 
is reproduced, whereas xq becomes one-loop order for 
6 < d < 8, see ([M]) . 

Below 6 dimensions, the leading scaling behaviour can 
be obtained by using fixed point values in ([8]), and also 
neglecting correction terms providing: 

w*h = \{H/kTf - {m\ - t) iw*q) + ^{n - 2) iw*q)^ 
mi ^ m\ — T + [w* q), 7712 = —(w*g), 7713 = 0, 
wi ^ w* , Wi = Q, 7 = 2 ... 8. 

After eliminating r by the equation of state, we are left 
with a two-parameter theory, with the simple RG flows 
close to the fixed point: 

Q~ (2- e/2 + 7772)9 and (H/kT)^ Xa (H/kT)^, 

(16) 

with Ao = 4 - e/2 - rj* /2 and ?y* = -e/3. The scaling 
fields can now be expressed as 

g., ~ {w*qy' g,{x), 

(H/kTf , . 4-e/2-77*/2 
X = ^ ' , ; and d = ' ' ' 



FIG. 1: AT lines on the two sides of the upper critical dimen- 
sion 6. See Uni, (HHJ and footnote [g] 
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(a) d > 6 



(b) d < 6 



From (|16p follows that x is invariant under renormaliza- 
tion, and Zi must be (2 — e/2-|-77*/2)"^ for gi = Xigi be 
satisfied. Any observable O satisfying the approximate 
RG flows O ~ fco O around the fixed point can now be 
written as O (y;*q)fco/(2-e/2+J7'/2) times a function of 
X. For a mass fc = 2 - 7;* ((5 - 1) (2 - e/2 -I- 7y*/2), and 
therefore the replicon mass takes the scaling form 



(u.*g)*-ifR(x), X 



{H/kTf 
{w*qY 



(17) 



The most important new feature of ([T7|) when compared 
with the d > 6 cases, Eqs. ([9]) and ([T3|), is the lack of the 
second scaling variable, which is proportional to 77. The 
AT line ends now in the zero field critical point ^ even 
for 77 small but nonzero: 

{H/kTf = xo {w*q)\ xo = -77 + [2 + 0(77)]e + 0{e^), 

(18) 

and it disappears completely for ti > 2e. 

To conclude, we followed the fate of the AT line from 
mean field down to d = 6 — e. An exceptional feature of 
the spin glass case (71 = 0) is that the runaway flows to- 
wards zero-temperature behavior — found below d = 6 
— originate in the close vicinity of the zero-field fixed 
point. Our results do not exclude a possible lack of the 
AT surface in d = 3 — as suggested by recent numeri- 
cal works 0, 01 — , a scenario, with some lower critical 
dimension to explain this, has been suggested in [l9| . 
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2 When 6 < d < 8 the AT line takes the form (H/kTf 
nl-2/£ _ Ti > 0, where qc ~ nr"^!^ . For d > 8 the 

mean field phase diagram restores 
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